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Recovery  of  Images  from  the 
AMOS  ELSI  Data  for  STS-33 


1.  INTRODUCTION 

During  (he  STS-33  space-shuttle  mission  in  kite  November  19$9,  n  number  of  instruments  of  the 
Air  Force  Maui  Optical  Station  (AMOS)  were  used  to  observe  five  overflights  of  the  shuttle  to  study 
the  interaction  of  spacecraft  contaminants  with  the  residual  atmosphere.  These  m^tsuiunenis  were 
part  of  a  secondary  shuttle  experiment  called  the  AMOS  Calibration  Tests. 

One  of  tire  instruments  used  was  lire  Enhanced  Longwave  Sjtectral  Imager  (ELSI),  which  is  an 
infrared-sensitive  array  detector  located  at  a  focal  plane  on  the  side  mount  of  the  1,6-mcicr  telescope. 
Data  from  this  instrument  are  recorded  on  tape  in  both  video  and  digital  formats.  Tire  ELSI  was  used 
on  three  passes,  orbits  21, 37,  and  67  on  24, 25,  and  27  November.  These  data  arc  normally  unclassified, 
bur  special  circumstances  required  the  ELSI  data  of  24  and  25  November  to  be  classified.  The 
unclassified  data  of  27  November,  in  video  format,  were  Juindcarried  to  Geophysics  Laboratory  (GL)  tit 
the  beginning  of  December  19S9;  the  classified  data,  in  digital  format,  were  sent  to  GL  from  AMOS 
some  time  in  February  1990. 

Since  the  Spacecraft  Interactions  Branch  (Pi  IK)  does  not  have  facilities  for  processing  classified 
data,  the  production  of  images  from  the  digital  tajies  was  undertaken  by  the  Data  Systems  Branch 
(LCY).  PUR  provided  a  written  description  of  the  taj>e  format,  based  on  information  provided  by 
AMOS,  and  LCY  developed  the  routines  needed  to  read  the  tajies,  but  when  the  data  failed  to  yield 
recognizable  images,  the  effort  was  at  an  impasse.  To  j>crmir  PHK  to  look  Into  the  problem,  it  was 
decided  to  work  on  the  unclassified  data  of  orbit  67,  which  carries  the  test  number  ATN  47434.  We 
requested  and  received  (on  23  March  1990)  a  one-minute  sample  of  digital  data,  corresponding  to  the 
interval  during  which  the  best  images  were  seen  on  the  videotape. 

Examination  of  these  data  showed  that  the  instrument  was  suffering  from  several  malfunctions, 
one  of  them  serious  enough  to  make  the  data  unusable  when  processed  in  the  usual  manner.  This 
report  describes  the  problem  and  the  development  of  a  method  to  deal  with  it,  specifically,  the 
construction  of  an  algorithm  capable  of  effecting  a  complete  recovery  of  the  images. 

(Received  for  publication  19  April  1990) 
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DESCRIPTION  OF  THE  INSTRUMENT 


2, 


The  EtSI  instrument  employs  two  array  detectors;  a  mid-wave  (MWIR)  detector,  which  is  a 
gallium-doped  silicon  array  of  45  by  90  pixels,  and  a  long-wave  (LW1R)  detector,  whieh  is  a  HgCdTc 
array  of  64  by  64  pixels,  T).c  image  can  be  moved,  by  a  switch  in  the  control  room,  from  one  detector  to 
the  other.  The  spectral  range  and  larger  field  of  view  of  the  LWJR  detector  make  it  more  useful  for 
acquiring  a  target;  after  acquisition  the  image  can  be  switched  to  the  MWIR  detector.  The  recording  for 
orbit  67  includes  data  from  both  detectors,  as  noted  below,  but  only  the  MWIR  data  arc  of  interest  here. 

The  ELS l  uses  a  rapid-ehopping  scheme  to  remove  the  instrument  and  sky  background  against 
which  the  traeked  object  is  seen.  The  image  is  recorded  first  with  the  object  in  the  chosen  location  on 
the  array  detector  and  then  again  after  *‘je  object  has  been  deflected  slightly  by  moving  one  of  the  flat 
mirrors  in  the  optical  path.  Subtract  the  second  image  from  the  first  is  extremely  effective  in 
producing  a  elcar  image  of  the  traeked  object  against  an  almost  black  background.  The  object  is  not 
deflected  enough  to  move  it  completely  off  the  detector,  so  the  subtraction  results  in  a  positive  image  at 
the  chosen  location  and  a  negative  image  at  the  deflected  position;  the  latter  is  removed  by  setting 
negative  pixel  values  to  zero,  leaving  only  a  patch  of  completely  black  background.  The  dcfleetion, 
which  is  called  nodding,  is  done  at  a  constant  nominal  rate  of  40  cycles  per  second,  implying  an  image 
rate  of  SO  per  second. 

Analog  data  from  the  detector  arc  digitized  by  equipment  mounted  on  the  telcseope  and  then 
transmitted  to  the  control  room,  where  they  arc  recorded  on  digital  tape  ‘Ihe  recording  process  is 
further  described  in  the  diseussion  that  follows.  The  digitized  data  are  also  sampled  and  reconverted  to 
analog  form,  in  a  standard  video  format,  for  display  on  a  video  monitor  and  recording  on  videotape, 

3.  TATE  FORMAT  AND  CONTENT 

ELSI  digital  data  arc  on  standard  9-lrack  6250-hpi  tapes.  Details  of  the  tape  format  arc  shown  in 
Figure  1.  Physical  records  are  unformatted  binary  and  contain  16  kilobytes  (16,384  bytes),  except  for  the 
>st  record  in  a  file,  which  may  be  short.  A  tape  may  have  one  or  more  files,  terminated  by  cnd-of-file 
marks,  representing  different  time  intervals.  Each  file  begins  with  a  file  label  of  1024  bytes,  followed  by 
an  integral  number  of  ELSI  logical  records  of  185,344  bytes  each.  Therefore,  there  arc  more  than  1 1 
physical  records  for  each  logical  record,  and  boundaries  of  physical  and  logical  records  are  not  simply 
related.  Each  logical  record  contains,  in  this  order,  an  8-byte  spacer,  twenty  images,  and  a  1016-bytc 
record  libel. 

The  present  purpose  requires  no  information  from  the  file  label  or  the  8-byte  spacers.  The  only 
item  needed  from  the  re  *>rd  labels  is  the  time,  which  is  given  in  units  of  seconds  and  milliseconds  after 
the  reference  time  of 0000  GMT  on  01  January  1970.  The  seconds  component  is  an  unsigned  four-byte 
integer  in  bytes  21-24;  the  milliseconds  component  is  an  unsigned  two-byte  integer  in  bytes  25-26.  Pixel 
values  in  the  image  arrays  are  unsigned  two-byte  integers.  Integers  are  written  to  tape  starting  with  the 
most  significant  byte. 

Each  ELSI  image  contains  9216  bytes  or  4608  pixel  values  in  a  6x16x48  image  array.  These  arc 
relitcd  to  the  45x90  MWIR  detector  array  as  illustrated  in  Figures  2  and  3.  Figure 2  shows  how  the 
detector  array  is  laid  out  in  six  scctior*  that  have  separate  readouts.  Data  are  recorded  on  tape  in  the 
order  in  which  they  come  off  the  detector.  The  readout  takes  a  value  from  the  same  location  on  each  of 
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Figure  I.  Tape  format  used  for  the  recording  of  ELS'I  data 


the  six  sections  before  it  proceeds  to  the  next  location  in  the  next  row.  It  progresses  through  all  row 
locations  in  a  column  before  moving  to  the  next  column.  Thus,  the  resulting  image  array  has  dimensions 
of  section,  row,  and  column  (s,r,c),  where  the  s  index  cycJcs  most  rapidly  and  the  c  index  most  slowly. 
There  is  a  further  complication  in  that  the  readout  progresses  through  16  (not  15)  rows  in  each  section 
and  48  (not  45)  columns.  The  resulting  spurious  values  are  referred  to  as  overclocked  pixels  and  must  be 
discarded  in  processing  the  data.  These  are  apparent  in  Figure  3,  which  illustrates  the  structure  of  the 
image  array,  the  ovcrclocked  pixels  are  represented  by  open  squares,  while  real  pixels  arc  shown  as  solid 
squares. 
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4. 


EQUIPMENT  AND  PROCEDURE 


Figure  4  show*,  m  a  function  of  time,  the  data  recorded  for  orbit  67,  Culmination,  the  time  of 
closest  approach  of  the  shuttle  to  AMOS,  k  Indicated  by  an  arrow  at  the  top.  Two  tracks  of  solid  bars  at 
the  bottom  show  which  detector  was  In  use  at  a  particular  time.  Above  that  is  shown  the  one-minute 
period,  just  after  culmination,  for  which  digital  data  were  requested;  It  ^eludes  a  small  amount  of  L.WIR 
data.  The  data  actually  received  arc  Indicated  at  the  top;  they  cover  the  requested  period  and  continue 
until  the  tape  reel  is  filled.  The  division  Into  two  files  probably  reflects  a  switchover  to  a  new  tape  unit 
during  the  original  recording.  It  provided  a  cottvenknt  starting  point  for  our  examination;  all  of  the  data 
discussed  below  come  from  file  2. 
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Figure  4.  Summary  log  of  the  orblt-67  ELSI  data  used  in  this  work 


The  equipment  used  was  a  Zenith  Z-248  desktop  system,  a  PC-AT  clone,  equipped  with  dual  22- 
megabyte  Bernoulli  drives  in  addition  to  20-  and  40-mcgabyte  internal  disks.  Also  included  was  a 
Cordata  LP30QX  laser  printer  equipped  with  a  custom  font  of  12S  grey-scale  characters.  Except  for  fad: 
of  a  tape  drive,  this  system  proved  nearly  Weal  for  examining  these  data.  Data  were  transported  to  this 
system  on  a  Bernoulli  disk.  The  first  six  megabytes  of  tape  file  2  were  copied  from  tape  to  a  Microvax 
user  dirccioiy  on  disk  and  then,  via  local  network,  to  the  Bernoulli  cartrWge.  This  data  extract 
comprises  a  little  more  than  32  ELSI  records.  Figure  5  indicates  schematically  how  the  data  treatment 
progressed  from  this  point  through  subsequent  steps  described  below. 
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Figure  5.  Schematic  diagram  of  the  processing  procedure 


The  approach  to  the  problem  set  four  sequential  objectives:  (I)  to  discover  what  was  wrong  with 
the  ELS  I  data  or  the  previous  efforts  to  use  it;  (2)  to  see  whether  the  data  of  orbit  67  could  be  processed 
to  produce  images  that  matched  the  video  data;  (3)  to  transfer  whatever  method  was  developed  to  the 
LCY  facility,  applying  it  first  to  the  orbit-67  data  for  confirmation;  and  (4)  to  use  the  established 
technique  to  process  the  classified  data. 
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EXAMINING  THE  DATA 


Pixel  values  in  (he  Et£l  daw  originate  in  a  12-bit  A/D  converter  and  are  pmlttvc  integers  ranging 
from  «ro  through  4095;  this  is  true  for  all  twenty  of  the  raw  images  of  an  ELS1  record.  Values  arc 
stared  as  two-byte  integers  with  the  four  most  significant  bits  set  to  zero.  In  practice,  the  values  may  be 
confined  mainly  to  a  much  narrower  range  than  the  full  4096  values  allowed.  The  principal  problem 
with  the  data  was  revealed  immediately  when  the  distribution  of  values  in  a  small  sample  was  examined. 
All  values  lay  between  3456  and  35S3  or  between  3966  and  4095.  What  these  ranges  have  in  common  is 
that  bits  11, 10,  $,  and  7  of  the  16-bit  integer  always  have  the  value  one;  i.c.  these  bits  are  'stuck"  in  the 
on  state.  Bits  are  numbered  from  tSB  «0  through  MSB*  15,  A  closer  look  abo  showed  that  bit  I  was 
also  Stuck  in  the  m  state.  With  five  of  twelve  hits  inoperable,  the  data  arc  enormously  ambiguous.  The 
resulting  response  curve,  plotted  in  Figure  6,  shows  that  each  recorded  number  represents  wh  the  one 
true  value  but  an)*  one  of  sixteen  possibilities.  The  scale  of  this  curve  docs  not  show  the  effoct  of  stuck 
bit  #1,  which  has  a  data  value  of  only  2.  As  will  be  seen  later,  its  effeet  on  an  image  will  be  negligibly 
small,  so  it  is  ignored  in  alt  of  the  following  discussion. 


Figure  6.  Plot  of  recorded  datum  value  versus  the  true  pixel  value 
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This  discovery  was  reported  to  AMOS  with  a  request  to  cheek  the  current  status  of  the 
instrument.  We  were  subsequently  informed  that  the  bits  were  found  to  be  stuck  exact!)*  as  reported. 
We  were  not  told  the  location  and  cause  of  the  malfunction,  but  because  the  video  data,  which  tire  stated 
to  be  derived  by  reconversion  from  the  digital,  seem  to  be  unaffected  by  the  stuck  bits,  we  suspect  lire 
bits  are  stuck  in  the  recording  equipment,  not  in  the  data  from  the  telescope. 

The  next  step  was  to  write  routines  to  rccogntee,  read,  and  manipulate  ELSI  records,  as  dictated 
by  the  format  of  Figure  1.  One  routine  was  required  to  clean  out  extraneous  material,  because  in 
rending  the  tape  the  Micros?.*  copy  utility  hod  inserted  an  cnd-of-record  mark.  ODOA,  after  the  data 
from  each  tape  record.  The  next  routine  operates  on  each  ELSI  record  to  decode  time,  di.Mr.trd 
ovcrclockcd  pixels,  reorder  real  pixels,  and  write  a  new  data  file. 

Examining  the  ditto  with  these  tools  disclosed  two  more  apparent  faults.  First,  one  sect  kin  of  il>c 
detector,  fifth  from  the  top,  was  completely  dead,  with  all  pixels  reading  the  saturation  value  4095.  This 
is  a  relatively  minor  problem  in  this  case,  since  the  tracked  object  was  kept  well  within  the  upper  lour 
sections  of  the  detector.  Second,  time  values  recorded  in  successive  frame*  do  not  advance  hv  the 
nominal  amount  of  250  ms  (20  frame*  at  a  rate  of  80  per  second).  Figure  7  shown  the  stdvanec  of 
recorded  time  over  a  series  of  .12  frames,  revealing  that  the  time  interval  per  frame  oscillate*  irregularly 
between  about  300  and  700  ms,  with  the  average  being  about  *410  ms.  It  is  mil  known  witeiher  this 
variation  is  real  or  the  result  of  inaccurate  reading  of  tltc  dock.  It  too  is  a  minor  ptohlem  here,  since 
accurate  timing  is  m>t  required.  These  faults  have  not  yet  been  reported  to  or  confirmed  by  AMOS 
personnel. 


3  *  i  S  T  t  S  IO  II  1313  H  l»  I*  IT  l«  I!*  23  31  »53  S«  ifi  ST  C«  S3  30  31  M 
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Figure  7.  Observations  of  recorded  time  versus  the  record  number 
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A  general-purpose  image- processing  routine  wu  written  to  permit  the  production  of  hardcopy 
grey-scale  prints  from  the  data;  this  was  needed  if  any  method  of  salvaging  the  data  were  to  be  found.  In 
its  latest  form,  the  routine  incorporates  a  number  of  sequential  but  repeatable  steps,  as  follows: 

a.  /fend  the  nest  record,  interchange  bytes,  and  convert  to  integers. 

b.  Subtract  each  odd  background  image  from  the  following  even  image 

c.  Correct  the  data  using  an  incorporated  algorithm. 

d.  Average  over  a  specified  number  of  successive  image  pairs. 

c.  Process  the  composite  image,  even-minus -odd  or  odd-minus-even,  by  finding  maximum 
and  minimum  data  values,  setting  and  adjusting  threshold  and  saturation  levels  for  the 
grey-scale,  and  converting  pixel  values  to  grey-scale  values. 

f.  Output  the  resulting  image  to  a  print  file  for  transfer  to  the  laser  printer. 

A  bateh  file  was  used  to  copy  the  prim  file  to  the  laser  printer  along  with  the  required  grey-scale 
setup  and  reset  commands.  With  defect-free  data,  these  procedures  wilt  produce  a  good  image.  The 
custom  GREY3  font  relates  coverage,  or  density,  of  a  square  pixel  to  an  ASCII  character  number  by  a 
conversion  curve  that  consists  of  three  connecting  linear  segments  that  lie  about  midway  between  a  pure 
linear  and  a  pure  logarithmic  scale. 

6.  CHOOSING  A  CORRECTION  METHOD 

A  first  version  of  the  image  processor,  with  no  correction  algorithm,  was  applied  to  the  unaltered 
data  to  see  how  much,  if  any,  of  the  true  image  could  be  seen  despite  the  stuck-bk  malfunction.  Figure  $ 
shows  the  result  of  applying  this  normal  processing  to  the  first  image  pair  of  the  first  record.  The  data 
printed  represent  image  2  minus  image  1,  with  negative  pixel  values  set  to  zero.  The  grey-scale 
threshold  (black)  was  set  at  zero,  and  the  saturation  (white)  was  set  to  the  maximum  value  found  in  the 
data. 

This  image  reveals  several  interesting  features.  First,  the  fact  that  detector  section  5  is  not 
operating  is  evidenced  by  the  black  band  across  the  lower  half  of  the  image.  Second,  the  location  of  the 
true  image  of  the  space  shuttle  fa  easily  discernible  as  the  light  area  centered  near  the  left  edge  of  section 
2,  roughly  at  column  10,  row  27.  The  light-to-medium  grey  pixels  clearly  show  the  shuttle  outline  as  it 
appears  in  the  videotape.  The  video  image  of  the  shuttle  b  a  vaguely  triangular  blurry  white  silhouette 
against  an  almost-btack  background;  there  b  some,  but  not  much,  variation  within  the  silhouette.  In  the 
digital  image  roughly  half  of  the  pixels  within  the  silhouette  have  dropped  out  and  appear  black.  At  the 
location  to  which  the  image  b  deflected,  centered  roughly  at  column  10,  row  45,  there  b  a  collection  of 
bright  pixels  forming  a  sort  of  ghost  image  where  there  should  be  only  a  completely  black  area  or  “hole". 
A  number  of  other  false  bright  pixels  are  scattered  about  the  rest  of  the  background. 


(on  facing  page) 

Figure  8.  A  background-subtracted  image  produced  by  normal  processing 
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Figure  9  shows  image  2  alone,  without  subtraction  of  the  background  image.  There  is  no 
recognizable  structure  in  the  region  of  the  shuttle  silhouette.  The  pixel  values  show  some  cohesive 
appearance;  from  upper  left  to  lower  right  the  pixels  successively  assume  tones  of  dark,  light,  and 
mediant  grey  in  three  rather  distinct  areas.  The  presence  of  only  two  boundaries  suggests  that  no  more 
titan  two  of  the  larger  jumps  in  the  response  curve,  3S4  or  6-10,  are  encountered,  and  only  one  if  the 
second  boundary  Is  a  retracing  of  the  first.  We  therefore  expect  the  distribution  of  true  data  values  to  be 
confined  to  a  very  narrow  part  of  thcv.cro-to-1095  range. 

In  the  background-subtracted  image  of  Figure  8,  the  abrupt  disappearance  of  grey  tones  or  the 
appearance  of  bright  spots,  as  one  moves  spatially  across  the  image  along  some  line,  is  easily  understood 
on  the  basis  of  the  response  curst  of  Figure  6,  It  occurs  when  the  true  value  of  the  even  (odd)  image 
crosses  one  of  the  discontinuities  in  the  response  curve  while  the  odd  (even)  does  not.  The  dark  pixels 
within  the  shuttle  outline  are  caused  by  a  negative  jump  in  the  even  image  or  a  positive  jump  in  the  odd 
image.  Hie  bright  pixels  within  the  blaek  hole  are  caused  by  the  reverse  condition.  There  arc  relatively 
few  fake  bright  pixels,  and  presumably  a  similar  number  of  false  blaek  pixels,  throughout  the 
background  because  it  is  rather  unlikely  that  the  illumination  on  such  pixels  will  he  much  different  in  the 
even  and  odd  images  and  therefore  unlikely  to  straddle  a  jump  in  the  response  curve. 

Tire  fact  tlrat  many  of  the  false  pixel  values  cm  be  recognised  by  visual  examination  suggests  that 
an  algorithm  aright  be  devised  to  correct  the  data  despite  the  ugly  nature  of  the  response  eurve.  Tire  eye 
discerns  bad  pixels  because  of  their  sharp  jumps  from  the  values  of  neighboring  pixels,  lire  jumps  arc 
noticeable  because  of  tire  relatively  low  spatial  gradient  of  intensity  within  the  true  inrage.  In 
approaching  the  shuttle  silhouette  from  any  side,  we  notice  that  the  intensity  increases  smoothly  over 
several  pixels  before  hitting  a  jump  tlrat  makes  it  false.  This  suggests  tlrat  the  average  spatial  gradient 
may  be  much  less  than  I2S  units  per  pixel-spacing  distance.  Since  the  smallest  jump  in  the  resporrse 
curve  is  I2S  units,  these  jumps  may  well  be  recognizable  in  goiirg  from  pixel  to  pixel  in  the  image.  As 
noted  above,  tire  true  pixel  values  seem  to  be  confined  to  a  narrow  portion  of  the  full  12-bit  range. 

lire  foregoing  observations  suggest  using  a  very  simple  and  straightforward  algorithm  for 
corrcetrng  the  data.  This  method  takes  one  starting  pixel  as  a  reference  point  and  then  moves  along 
multiple  tracks  across  the  detector,  adjusting  the  value  of  each  pixel  according  to  whether  and  how  much 
u  Iras  jumped  relative  to  the  preceding  pixel  on  the  track.  Tire  pattern  of  tracks  is  a  tree  with  enough 
branches  to  cover  every  pixel.  Since  any  error  in  identifying  a  jump  will  affect  every  subsequent  pixel  on 
a  branch,  it  is  well  to  make  the  branches  ns  short  as  possible.  It  is  not  really  necessary  for  the  starting 
pixel  to  have  a  true  value,  since  the  absolute  level  of  the  result  Is  not  important  unless  successive  images 
are  to  be  averaged.  But  in  any  ease,  the  amount  of  its  error  will  be  apparent  from  the  difference  in  the 
resulting  distributions  of  positive  and  negative  values;  the  black  background  will  not  be  roughly  zero  but 
some  value  that  approximates  tire  $iz.c  of  jump  by  which  the  starting  pixel  was  in  error,  and  the 
magnitudes  of  the  most  positive  and  most  negative  pixels  will  differ  by  twice  the  size  of  that  jump. 


(on  facing  page) 

Figure  9.  A  single  even-numbered  image  produced  by  normal  processing 
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7. 


CONSTRUCTING  THE  ALGORITHM 


Far  clarity  In  the  following  discussion,  all  results  presented  will  be  for  the  Image  pair  already 
shown  (file  2,  record  1,  image  pair  1),  and  references  to  the  sample  image  will  mean  the  background- 
subtracted  image  unless  stated  otherwise.  Also,  in  mueh  of  the  work  to  follow,  only  the  upper  two- 
thirds,  rows  l  through  60,  of  the  detector  will  be  used.  The  fifth  section  is  dead,  and  since  the  detached 
sixth  section  is  of  little  use  by  itself,  it  was  usually  set  to  zero. 

We  know  that  the  distribution  of  values  in  a  single  image  is  confined  to  the  two  ranges  noted 
earlier  3156-3583  and  3%$-4095,  but  this  will  be  different  in  a  background-subtracted  image.  Figure  10 
shows  the  distribution  for  the  sample  case.  This  and  the  following  figure  use  logarithmic  scales  except 
that  zero  has  been  placed  at  the  0.1  level,  Values  are  less  eonfined  because  two  quantities  are  involved 
in  the  subtraction.  We  may  reasonably  guess  that  many  of  the  values  below  12S  are  true,  while  most  of 
the  two  higher  grottps  arc  false.  We  also  note  that  over  most  of  the  image,  where  there  is  no  tracked 
object,  both  the  even-  and  odd-image  values  of  a  pixel  may  be  grossly  affected  by  the  sluck-bit 
malfunction  while  their  difference  remains  correct.) 


MAGNITUDE  OF  VALUE 


Figure  10.  Distribution  of  pixel  values  in  an  uncorrccted  image 
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The  distribution  of  the  nctu.il  pixel-to-pixcl  jumps  In  «n  image  is  most  important  to  know.  Figure 
11  shews  the  distribution  of  jump  sizes  encountered  in  scanning  left  to  right  along  all  of  the  first  60  rows 
of  an  image.  Jumps  having  magnitudes  of  l,  3,  and  5,  in  units  of  128,  arc  expected,  since  they  occur 
directly  in  the  response  curve;  the  pixel  needs  to  jump  in  only  one  image  of  the  subtracted  pair.  We  also 
sec  many  jumps  of  magnitude  4.  In  this  ease  the  pixel  must  jump  In  both  Images,  for  example  { +  3)  ( 1) 
or  ( *  5)-(  +  1).  Both  of  these  combinations  arc  likely,  as  cart  be  seen  by  juxtaposing  the  curve  of  Figure 
6  and  its  inverse;  a  difference  of  only  12S  in  the  true  values  of  even  and  odd  images  can  cause  these  jump 
combinations  to  occur  if  the  spatial  gradient  has  the  same  direction  in  both  images,  as  is  usually  the  case. 
Conversely,  jumps  of  magnitude  2  and  6  are  almost  never  seen.  These  also  require  jumps  in  both 
images;  for  example  ( +  3)-(  + 1)  or  ( +  5)-(-l).  Sueh  jumps  are  also  possible  with  the  same  even-odd 
difference  of  12S,  but  the  spatial  gradients  in  the  even  and  odd  images  must  be  in  opposite  directions. 
Since  this  is  very  unlikely,  these  jumps  are  very  rare. 


Q 


s 


SIZE  OF  JUMP 


Figure  11.  Distribution  of  pixel-to-pixel  jumps  in  an  unconnected  image 
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Clearly,  the  important  feature  ef  this  distribution  tlta<  makes  correction  possible  is  the  clear 
separation  of  the  magnitude-0  and  magnitude*!  groups.  Tills  repeals  that  the  pixcl-to-pixcl  variation  in 
the  true  image  does  not  exceed  about  <i§,  white  a  s<uek*bit  jump  os-era  diseontlnuity  of  12S  never  results 
in  a  jump  of  less  than  6-t.  Similarly,  jumps  of  magnitude  3,  *1,  or  5  arc  equally  unambiguous.  This  implies 
that  an  algorithm  based  on  jump  size  alone  will  suffice,  and  there  will  be  no  need  to  invoke  such 
additional  criteria  as  the  local  trend  of  the  spatial  gradient,  cte. 

The  correction  algorithm,  therefore,  was  constructed  to  work  as  follows.  The  starting  location  is 
taken  to  be  the  top  ef  the  middle  column  (eolumn  23,  row  i)  and  the  first  track  runs  down  column  23  to 
row  60.  From  the  central  column  ef  corrected  pixels,  each  subsequent  track  runs  along  a  row  from  the 
central  eolumn  to  the  right  or  left  edge,  so  the  complete  pattern  comprises  12!  tracks,  with  only  the 
initial  track  covering  more  than  23  locations,  At  each  step  along  the  track  the  observed  jump  to  the 
current  pixel  is  classified  as  to  the  discontinuity  included,  an  integral  number  between  *6  and  4-6.  Tills 
number  Is  added  to  the  cumulative  total  from  the  beginning  of  the  track  to  find  the  amount  by  which  the 
current  pixel  should  be  adjusted.  The  cumulative  adjustment  for  each  pixel  Is  stored  in  a  separate  army, 
and  the  pixels  are  adjusted  after  completing  all  iraeks. 

The  threshold  levels  used  to  determine  jump  size,  that  Is,  the  location  of  boundaries  between  the 
groups  in  Figure  !!,  were  allowed  to  be  individually  adjusted,  but  in  practice  setting  them  all  to  the 
appropriate  multiple  of  (hi  was  found  to  be  adequate. 

It  may  be  useful  to  note  that  we  first  tried  to  use  a  similar  algorithm  applied  separately  to  the  odd 
and  even  images  of  a  pair.  It  was  thought  that  this  might  prove  best,  so  that  legitimate  variations  in  both 
images  would  not  be  added,  possibly  making  their  sum  large  enough  to  be  eonfused  with  a  magnitude*  l 
jump.  Tills  was  a  mistake.  That  effect  was  outweighed  by  effeei  of  inherent  pixel-to-pixci  variations  in 
response.  The  latter,  which  creates  the  need  for  flat-fielding  techniques  in  the  first  plaee,  proved  to  be 
larger  than  pixel* to-pixe!  variations  in  the  image.  They  were,  in  faet,  so  large  that  they  could  not  be 
distinguished  from  magnitude*  I  jumps,  and  the  algorithm  failed. 

8.  DESCRIPTION  OF  RESULTS 

The  algorithm  proved  completely  successful  in  correcting  the  orbit-67  data  used  in  its 
development.  Figure  12  shows  the  fully  corrected  sample  image.  In  this  ease,  grey-scale  saturation 
(white)  was  set  to  correspond  with  the  brightest  pixel,  and  grey-seale  threshold  (black)  was  set  somewhat 
above  zero,  as  is  usually  done,  to  suppress  the  empty-sky  background  in  addition  to  eliminating  the 
deflected  image  due  to  nodding.  Figure  13  sho**  N:  same  image  without  the  sky  background 
suppressed  (grcy-seale  threshold  m  zero)  '  Hack-ground  lias  a  fairly  uniform  mottled 

appearance,  and  the  black  hole  of  the  deflect  ’ble. 


(on  facing  page) 

Figure  12.  The  fully  corrected  image  with  sky  background  suppressed 
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Same  testing  of  the  results  was  danc,  but  no  significantly  different  data  were  available  far  a 
rigorous  evaluation.  As  long  ns  the  jump  distribution  shows  complete  separation  between  groups,  the 
correction  should  be  perfect.  U  is  not  an  approximation,  so  even  the  most  subtle  features  in  the  true 
image  will  be  preserved.  Murco-er,  if  an  error  in  detecting  and  classifying  a  jump  were  to  occur,  it  would 
expose  itself  by  producing  a  very  noticeable  defect  ft  ibc  resulting  image,  a  defect  propagating  along  the 
remainder  of  the  correction  track.  litis  was  demonstrated  in  testing  by  deliberately  forcing  tin  error. 
\Vc  expect  tltat  such  a  defect  would  still  be  apparent  even  if  tlte  image  were  too  active  for  the  correction 
to  succeed. 

Tlte  successful  correction  of  these  data  reveals  one  more  condition  that  may  be  another 
significant  fault  in  the  instrument.  The  trucked  objeet  seen  in  Figures  12  and  13  is  the  space  shuttle, 
which  is  only  -10  meters  la  length,  at  a  distance  from  the  teleeope  of  more  than  200  miles.  Still,  in  Figure 
13  we  can  clearly  see  that  the  black  hole  of  tlte  deficetcd  image  has  obscured  a  small  pari  of  the  desired 
image.  If  tlte  object  were  either  larger  or  closer,  more  of  the  image  would  he  obscured.  Since  the 
background  subtraction  works  very  well  with  the  current  amount  of  image  deflection,  it  should  be 
po'Mblc  and  useful  to  increase  the  deflection,  moving  the  black  hole  farther  out  of  the  way.  We  have 
not  vet  enquired  about  the  practicality  of  doing  so. 

At  the  time  c»f  this  writing,  the  algorithm  has  been  transferred  to  LCY  for  processing  of  the 
classified  data.  A  preliminary  test  on  their  system,  using  the  same  unclassified  data,  appears  to  have 
produced  tlte  same  result,  but  because  of  security-classification  dilficultics  that  system  currently  lacks  -t 
grey -scale  output,  a  problem  that  is  currently  being  addressed. 

FORTRAN  source  listings  for  most  of  the  software  written  Tor  this  work  are  reproduced  ia  the 
A|  |  Ciuliv.  Only  the  most  recent  versions  arc  included.  Some  of  the  results  presented  above  were 
produced  with  earlier  versions. 


(on  facing  page) 

Figure  13.  The  fully  corrected  image  with  grey-scale  threshold  at  zero 
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Appendix 


software  source  listings 

FORTRAN  source  code  for  the  following  routines  ore  listed  below,  in  the  order  shown  here. 
The  eede  is  written  for  an  IBM  PC-eomp.liibIc  computer  running  under  DOS. 

DD  A  simple  dump  of  the  data  Hie,  in  512  byte  bteeks,  to  tire  monitor  or  a  disk  file. 

RR  A  routine  to  femove  spurious  end-of-reeerd  marks  inserted  after  each  tape 
record  by  the  tape-copy  utility. 

BE  An  ELS  I -record  reader  to  extract  and  examine  one  record  (spacer,  images,  and 

label)  at  a  time. 

TF  A**  intage  translator  to  read  the  raw-data  file,  extract  and  deeode  time,  discard 
ovcreloeked  pixels,  reorder  the  array,  and  write  a  new  data  file, 

DS  A  routine  to  compute  distributions  of  pixel  values  and  pixel-to-pixel  jumps. 

GG  An  image  processor  to  read,  subtract,  eorrect,  average,  process,  and  output  the 

results  from  one  ELSl  record. 
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RFAOlWJC 

lF(lsf.(i£.6)G0IQ  « 

.. — — secede  tie* 

»* 

IT»«60(!T.6C) 

KT  J).|T 
KTMJ-1I2 

... - display  B-byte  spacer  and  Us* 

W1TE<*.B18)8  KT  r 
818  FCRJlAtf 

*  IX.'EIcHT-BYTE  SPACER:  *, 

*  IX.^tsf-^COfiO  TIME: 

♦  12,J:'fI2,*:*,12,',',IJ,  •  UT\/? 

- display  28  iwges  in  hex 

KR1TEK815)  s 
815  FCRJtATf 

♦  IX, 'Display hex  dusp  of  Images  (Y/N)?  IN]  \\) 
READ{»,7£3)Z 

IF(Z.n£.'Y  .AfO.Z.NE.'y'JGOTO  22 


24 


w  :<  im5*t.ja 
ui}t£<* «lt)l«A5 


j  Mntet  fw«c* '.J2.A 
X, ’Enter... 0:  Display  i 


!  Display  ihlt  l**4«!»A 
i  f*Ip  to  next  Iwm',/, 
i  Read  the  next  rtc# rs’. 


,v '  .^„£{  0,111  per*  *••• 

M'l*  •2,H8*!‘!!co>o  * 

tF(19CR?nj,ljCQIQ  >7 

gtifeig- 

fxl'f^cticre^cwn 
63  M  llKUwl.ti 


(b)  *A) 


CWTH4C 
. itvj  rtWS 

IX. ‘Snip 


!&!>' 


IX, 'Snip  CO}tKfl/N«Ji 

iWSSSSi cl 


usn/Next  iMje/Read  recerd/txll’. 


FfAM.’S’.tt.r.EO.’s'H'KM 
Fll.ro.'N'.M.Iio.'n'jGQto  :a 
Ffl.lO.’R'.W.i.fO.T*  COlO  48 
n|.t^'£'.«.Z.to.vicalD  55 

JJ—Cg  COllWtt 

fCSKAItlx/eejln  next  lM»e  , 

A> 

Ihz.fO.’S’.w.i.Eo.’s’iJwo 
i'li.tO.'r.W.Z.tO.'e'jMIO  55 
WKtiwe  , 

IfcJIS 

CONTINUE  * 

—output  a  KKJIfletj  ElSI  rtari 
not  leplt«nted 
— optional  exit 

IX 'Read  another  £151  record  (Y/N)?  (v)  •  ,\) 
R|*0{*,7JCU 

IfU.EO.'N*  .Ofl.Z.EQ.'n'JGQTO  42 
CONTINUE 
-"dose  files 
CONTINUE 

flle-wite  code  not  present 
GOTO  55 

—end  ol  file  or  read  error 
IFflST.GT.BJCOTO  47 
WUEC.838 

FORKaTT  EXIT:  End-of-flle  encountered’) 

GOTO  42 

WITEKBJ2) 

FORHATt1  ASuRT:  Read  error  In  EISI.OAT') 

— nor o*l  exit 
CONTINUE 
DOSE  (50) 

UtlTE(*.899) 

FORMAT  >/. 


Read  error  In  EISI.OAT’) 


ENO  OF  EIS1  DATA  READER  (EE) 


TT  ROUTINE  (IMAGE TRANSLATOR): 


rsfcsw  IT 

—IMS  routine  reads  In  iMleal  reeards  f rm 
a  JUiWarfl  ttSi  data  S«d  W>HM  &ii  new 
logical  recem  that  remain  enly  the  if*«,  a 
reused  Identifier,  sod  UK  i-sf.ty  iMces.  The 
'l*«l  values  a^e  reordered  to  a  (wo*fli&eftsiena3 
array  of  <5  eeliws  by  $4  rpws  (free*  ins  tape 


i 


discarded  d-jrtna  reordering. 

,  i  ifftjm  i*  ii|mt  byui  aim. 

output  record  wfoui  h  I628ib  fcyies  (»iKi 
frisiM.sn  Ii6-b,te  header  and  M  tlfO-b>te 
ge*  Hi6  *  i iztih  »  16:316  ty;es  *  \l%y. 


The  input  record 
the  output  reenr 

peer?  - 


!*«•  variables,  etc  . 

pURAEURM  Atlg2t 
BiAStCiiSM 
CHASACllRM 

ousaeur.i 

OUSiCUS 
CHAStCU  R 
6iijfiGtlS*8 

im*  r»2 


|NIIG#H*4 

icy  vaifns 

tV  \AUN2 
W  v»i(v- 

vaum 
yau\~ 
s.  v*Usa 
..Cy  vaiins 
fey  viun: 
ievvsuNa 
-—define  constants 
(C.ya  2i  u?m  istf 


its 

rs»LsM>?*#  m* 

mic2i,R,fi«K 

litatr1 

If!  l!i5i}lji  ki!jSij!,|l 

ch  .wi.tc  i  ,«e ,  c  106) 
to  i  ,«{.  c  \  M  q  2  EJ 

MM- 

sal  * 


totm 


■  eie.ia.Moi 


tea 


•••display  weyaa  identifier 
(c-vjA  (}///////, 

i 


t  ZM 

•  tit 

«  let 

•  s?< 

•  w* 

•  2,’* 

•  MX 

•  M< 

•  mx 
•*  tax, 

•  mx 

i  m«; 

•  MX 

»  2i’< 

«  MX 

•  Mx 


ElSi  CATA-FJtE  IRWaATM  (IT) 

Re-order*  data  Into  on  array  ef 
IS  Cdjsrn*  fey  50  row*,  dlsturair 
cvereioviiea  pixels  in  rows  16  &r 
cash  section  and  eolusn*  <t*4R, 
Reads  fron  irput  Hie  US  j  .CAT 
Writes  to  output  file  Etsi.lKJ 
Cvth  flies  eust  already  exist 
In  the  default  directory. 

Press  return  to  begin  readlrj  ... 


3 


7  ca 


•  / 

.•v 

•7* 

•7’ 

•7’ 

•7’ 

»*/* 

•7* 

•7* 

•7 

•7 
■  •  «* 


602 


700)2 
FCATATiAl) 

—open  Irsut  and  output  files 
CPENiSD.FUE^ElSI  .DAT*  .FfAH.'BHiARY1 , 

access.  1 'secicnt  i  al  * ,  leaf  ai.  ist. status- 'cto1 ) 

Cy£N(tO.FIlE.*ELSI  .IKS' ,FCfi«.,fllKARr , 

ACCESS.'SEtjyENUAl’.liSIAT.JST.STATUS.'OtO') 

—read  and  display  file  label 

READiSOJA  ^ 

tKlTtt*,802)A 

FCf^ATt 

SOX.'Iaue  file  label  reads  ...*,/, 
T6(§X.64Al./J. 

20X,'Press  return  to  continue  ...') 
REAfi{*.700)Z 


•process  52  ELSI  records  (nt/sber  In  file) 

mi?E(»  ee?)65^'10^  byUs)/,855<<  *  52-« 
804  F0RKAT(‘8l 

■»  20X, ’Pause  before  each  record  (Y/N)7  [NJ  '  ,\) 
R£Ao(«.700)Z 
UUTE(".805) 

80?  FORKAT< IX,/) 

IPAUSE'B 

IF(Z.EQ. * Y* .OR.Z.EQ. ’y* ) IPAUSE-1 


26 


DO  38  l*£C«U2 
n 


M6  fow*l( 

*  aw, 

♦ 


GOTO  2 
C 


■  te  ff««j  fteard  M2, 
.rtlvfn  la  ecAiln* 


I)  IKCC 


Mi* 

*S)S!o3  ,n  M«rtMV 

’  .DGOJO  97 


mrmo” 

'Msfjse.Bjcoto  97 

Un« 

CO  J  1*1,4 

mfer1’ 


110 


M£ 

WJl*  I 

w 

••••aWplay  t!s< 

?&U'%lCXX 

*  JflX.'tlSI  ritccd'.U, 

♦  9;.!u»iY£fHi  TI«j  \ 

•  12*  i'.U.’J’.U.’.Mi,/) 

- MS  1M^‘I?l3y  **** 

s  11  ii»eiw«  * 

H  1  •  'flSI 
H  2  *  *nJiR 

11  j :  »• 

II  l :  3j8  !S 

w 

S*Dl*.*X 

see  mnar1 

msur 

hh1;.’K\* 

hoi)**  KiKi£R  • 


JIB 

314 

J15 

316 

318 


Hjljj-'OAIE  IS  • 
Hlj)*' 11/27/69' 

83ft  • 
ISM!!!.., 

FCRK*tf6x.  15) 

m 

wnE(ftj.ji«)i 

H  17j-'  PLUS  ' 
R*BIW* 

I*1MC*KT(4) 
URITE{R8,J18H 
FORf'Ald'l.lX,  ’KSEC') 


Hf 18J-R 
H  19>-SLAf*C 
H  28  ••0AV10  J.' 

H>  21  •'  KNECHT  • 

H  22  •'  3/26/98' 

«> 

t®l«8flX,to*r  rM> 

♦  17(16X6*8;/)/. 

♦  16X.48(lH-5,/) 


27 


lit  r 


ssgsftsr 


inj  recerdMV 


Iww 

- Ln,Wl*» 


guc^twlliwjjMooscr.uw.iea.imc) 

-wJU  n*w  rjjord  to  flit 

*11  ° 

•  out  retord'.IJ,' 

'i!  g|p)«wi./) 

. — djl  J*4"'*4 

115  [oimf(i5x!^H  record*  ereceited*) 

0010  tt 

* - —m)  of  fli*  or  re*‘ 

W«  ro«wf(i».trY'T'  ' - 

GOTO 

.8 


reed  error 

» 

‘EXIT*  EM  of  flit*) 


m  fowti ... _ 
cotort 

•5i  S«r*l(ilx!^EXlT:  Oitkflle  reed/wrlte  error') 
•••——rwool  wit 

Sr 


WlttJ* 

*19  foumf 

*  a*.:i 


ENO  OF  ELSl  OATA  TRANSLATOR 


./. 


3:3 


:>) 


DS  ROUTINE  (DISTRIBUTIONS): 

Mfindd?itrlbutlon  of  pixel  value*  and  Jiap* 
In  the  ten  eeeblned  (even-alnut-odd)  JUoes 
of  on  EISI  retold.  The  tile  of  the  dlttri- 
botton  bint  It  I. 

verilon  0000  -  0.  J.  Knecht  *1  April  11K 


VBflibltS 

CHARACTERS  AfU2ll6).0(U2M«KOS(l<W) 
CHARACTERS  K0r21,Z 
CHARACTERS*  SETH 
CHARACTER**  HUS2) 

E§uivALtNC£  |N^{ll!^ilj||(NisS}l]!oS(*l9J}) 
EQUIVALENCE  (0(l),«ll.l.I)J;(lK4.K6(U) 
-—open  filet 

0f£N(6d,FllE-'ELSI .  INS'  .FORN-'IINARY' , 

♦  ACCESS-  SEQUENTIAL', JQ$iAT*I0C,5TATU5*'01O') 
— eonttenti 
KMK-257J 
- Initialize 

isUBT-l 

ioisT-i 

IREC-I 

"affapr  * 


*  rqmiu///////. 

♦  2ix,- 


DISTRIBUTION  COWER  (08) 

Coaputei  the  dlttrlbutlon  of 
pixel  velutt  In  an  ELSI  laage 
Read*  input  froa  file  ElJTTmS 
lut  wist  on  default  drive) 
Write*  to  a  file  ipeclfiod  by 
keyboard  Input. 


:* 

•7' 

*7' 

•7’ 

•7' 

,'V 

•7' 

•7* 

if  o 
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.  .  SCUCt  OPERATION! 

R  »  jMd  an  Input  rtcerd 
5  *  5<*iraet  bkboreund* 
V  «  Values  dlitrlfutlon 
*'  *  Ju»ps  distribution 

t :  gjj{* 10  ■  «»* 


fsf 


sS£tE«10Ni- 
’atlen 


■>  W) 


''r 

.*y 

•7' 

•7* 

t»  t 


a» 

l.W.tt.Z. 

I.’SJ.OR,' 

l.'V'.OR, 

'.'J'.OR. 

g  .or. 

't'.OR. 


GOTO  4 
i  GOTO  11 
I  GOTO  jl 


I.M'  GOTO  _ 
!.  5*  GOTO  »1 
GOTO  W 


jl  record 
;?-<)G0!0  $ 


Ml 


It 

812 


wltEj*, 

f£WV.t{ 

'  cwrfwfc  rewd  m)7  CN3  'A> 

C£tT*»£CtAl0,N3,!*£C) 

ISU6T*a 

JOlST-d 
rOSTW 

2fX\ 

&a'.7M)Z 
GOTO  1 

-••subtract  within  each 
SFT|N^OTafO. ? JCOTO  12 

W I 


.'.[rm  return  to  continue 


(nm  pair 


12 


114 


41 

842 


FdRw|(iixtU  data  have  been  read.',/, 

.W  Si?"" 

GOTO  l 
ioisT-a 

,  mvm 

*  H5'!5ublr*ct,0°  ,s  «*pi«ud.'./, 

•  *E&  t"SS)J  “  •••'■> 

GOTO  1 

of  values 


1T|]S&  •  .DGOTO  42 
" — At{ 


Foamt{ 

MX, 'Subtractions  haw  nat  been  done. 

K*i(Jf7i*)2*lUrn  for  "#,n  ■«"“ 

GOTO  l' 


'.A 


.11  »!“> 

*  ooxjj' dlltrlbgtlon  {cr  ,n  P*lri  ’> 

moil!.!1'  ” 


ht»Mi.§ 

43  MSN  l  .| 

00  46  lPR.1,11 
K»1PR*2 
00  46  J-1,66 
00  46  1-1,45 
NV-hOfKJ.K) 

SSWT” 

ND5N(hNEG)-N0SN(*CG)*1 
GOTO  46 
NP0S-NV/I*1 


44 
46 

miTE(«. 

156  FOWAtf 

1  ?£K&lu*  distribution  completed. '.A 
+  Fare*™  t®  continue  ...\J 


afei*8 

GOTO  1 


61 


&&&****- 


862  Fomtc’ 

+  H5* '5obtr»ctloo*  have  not  been  done.  *,/, 
+  2 IX,  Press  return  for  aain  aenu  ...*) 


29 


Smckw*) 

ronwl(>, 

MXj'Coe^utJnj  distribution  for  *11  p*irs') 

fern’ . 


V)  44  |PR.),)I 

oo'w’j-i.ti 

W*G»W/I«1 
MSN(MCG)*OSN(»CG)tl 
GOTO  64 

PW{^)*C$P(NPCS)*1 

continue 

ttlTE(*,lfl) 

FORfiAtf 

JflX, 'Distribution  It  COBfiUtrt.',/. 
zex.'frttt  return  to  continue  ,»,',) 

—output  to  file 
IFf  lOlSI.tQ.UGOTO  12 

foumtf 

2*x, 'Distribution  Ms  not  been  computed. ', A 
2tx. 'frets  return  for  min  nenu  ...’) 

GOTO  I 

utllE(MM) 

rt«mt( 

28X, 'Printing  distribution  for 
MX. '*11  I«*§e  pairs  to  *  file.',/. 

MX. 'Provide  •  Tile  n*«* 

0PEN(77,FRE*'  ’.fORN-'FORNATlEO', 

STATUS^NEU' ,  ACCESS* 'SEQUENTI AL’ ) 
wjE|77.77lj(KBP(i).M,512).(MJSN(I).  1-1.512) 

• TEN  MAGES:  BINS  886  -  S12',52MH-),A 

■«liSl(S-,SiS,iS,:/jii?:  k  »: \i • 

li&ff  *  •  1X'75{,K*, 

OOSE  (77) 

GOTO  1  , 

—•nerael  exit 
CONTINUE 

mv?” 

ZBX,  *  - .  * ,/ 

2PX, ' I  ENO  OF  DISTRIBUTION  COMPUTATION •  / 

wb - J 

ENO 


SUBROUTINE  ROREC(A.D,NO,IREC) 

character*!  AfUiettho  uzitoj.oo.z 

INTEGER-2  Kl(45,96,26 
tAITE(-,382)IRt6 

382  F0RNAT(28X,  Reodlng  record*, 13,'  ...') 
REAO(U)A 

- switch  byte  order 

00  26  K-1.B1M 

!*(K-T)-2»1 

J-l-1 

OO-O(J) 

»il) 

28  CONTINUE 

”:-Kintd*t«  swple  (col  1-15.  r<w  1-fi') 

38(  F0RMAi('™ 

*  MX^'Print  dot.  tuple  to  file  (Y/N)7  (N)  '.\) 

188  F0WWT|A1)K 

IFfZ.Nt.'Y'.ANO.Z.NE.'y'lQOTO  48 
OPENtM.FilE-'  '  , FORM- 'FORMATTED' , 

♦  STATUS-’PCH') 


*  STATUS- ’f£H') 

DO  38  K-1,2 
00  25  J-1,68 

m  KKilt!i{S5,l“" 

25  CONTINUE 


.J.K). 1-1,15) 


m  »»! 


3  gas* 

DC 


oo  ?  ipa|r*i.i« 

00  r  ' 

00 

Si 

cm 

*Ui 
DO 


m 


BKaaaiwjtu 


GG  ROUTINE  (IMAGE  PROCESSOR): 


ffOOUfl  GC 

Recl^r5tf(w*l»plta*nttd  by  wrcitlrwite 
combined  mat.  Corrects  fiMti  ef  -4  tfireirth  «4 
U«i  ttf,  with  MM<rl(S  SVM  atcvt  and  6el{w 
these  v*lu<j,  MjtJ  m  distributions  obtained  with 
frOJTM  COO.  (0.  J.  Knecht  t)  April  lift) 


■CftritiH  plot  routine 
ef  it  reordered 
subtraction  and  t 
stucx  sit*: 
adjuii  aae,  aln,  — 


H 


.0(1)) 


k»w« 

l 

Mi 

♦ 

4 

♦ 

4 

♦ 

4 

♦ 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 


*«*$*  ln$n  CIS!  HU 
_  far 

. ™_ . -.v, 

MturiUon:  **  uriih  int  rtwlt  to  *  yey 
»t*}«  print  flit  fer  iht  1m*t  printer. 

— tkci^c  Mr  I  Mi  it . 

ewucTtft*) 

saw  »&!«•««.» 

VAita  _  . 

— open  flic* 

opihju.f  it£»'ctsi  :  ins*  lraw.*8JfWiv, . 

ACCESS* 'SCQUCNHAt,'IOSiAT*iCC,SIAT(£**CLD') 
—define  constants 

MUK-257J 

IB# 

-jjjjtj.il., 

ISWT*| 
lSAVt.1 
1KMM 
!«WC*4 
IREC-I 


5!ISi%KT 


TOW* 

MX,' 
MX  ' 
mx  • 
mx  • 
m  * 
mx  • 
2 «  * 
MX,' 
MX  ' 
MX, ' 
MX  * 
MX  ' 
MX  ' 
MX,' 
MX  ' 
MX,* 
MX,1 
MX  • 
MX  ' 
MX, ' 
MX' 


identifier 


Oman, 


ELS  I 


<«) 


mux  processor 
(Xevli I on  5.1) 
to  product  grey-scale  print! 
ftetdt  the  input  file  ElSUnG 
(•ust  wist  on  default  drive) 
writes  output  filet  naaed  in 
response  to  Interactive  query. 

SELECT  OPERATION: 

*  *  Read  an  input  record 
S  •  Subtract  backgrounds 
C  «  Correct  for  stuck  bits 
A  «  Average  selected  pairs 
P  •  Process  coapoiite  iaage 
0  ■  Output  to  print  file 


SELECTION:--)  *.\) 


\'ff' 

*■/' 

.'V 

.'V 

•V 

>y 

’7’ 

'7* 

•7* 

•7' 

if  « 
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r-*sl«t,»er»vlen 
!  «:$&« 


I 

,p  i 

« 


K&g. 

V£'.«. 


,'r* 

:4 


5310  4 
5310 
5310 

giS2! 

5310  it 


11 


.vjwto  i» 


t<u'Sj?WM.r.$.»i«TO  5 
« 


fowls' 

2*atx*0l*c<*rd  current  record  {V/NJ?  M  ’A) 


n 

112 

f 

12 


M-ii-.UO-it 

•iriasiidsawir ,r" 

data  have  been  read.' 
tct'frtit  return  fee  aaln  km  ...' 

Goto  1 

Oil  SUTfttW) 

151*1*1 


/. 


.14  »,,4) 

•  ^.'Subtraction  it  wepUita.*./, 

•  jei/fctM  rtturn  to  continue 

t«ltC<*,*22> 

152  feftttmcx.'fij  data  Have  betn  prepared.',/. 
«  J/x.'Pres*  return  fee  »4ln  *eno 

S  twi!;-,!'"'tc'""!' 

«  ttx, 'Correcting  data  of  rtcerd'.lJ,/, 

*  tax/lM  first  Jesse  it  rweeer  As,/) 
t>3  is  lFAlR»l,lfl  1 
KWIPAlft«l)*i*t 

K2‘K1*1 

Oil  KCtXINO.W.Kl) 

Oil  max  KJ.W.M) 
je  continue 
WE>N*S 
1M0C*» 

>»  r&m™ 

-  2£t, 'Correction  is  ccepleted.',/. 

’  Aii. 'Press  return  to  continue  ...',) 

ggHvwnz 

C010  1 

-9ver*3«  over  successive  pairs 
■U  JFStf#uf.EO.l.ASO,lSuar.t0\l)COIO  42 
taliE{».»42i 

842  FC8KiT(2SX,'No  data  have  been  prepared. 

♦  2«x.'?ress  return  for  eain  aenu 
RE^{;/(A1)-)Z 
C310  1 


./. 


■‘'•—specify  sussing 


42 

844  FCMIaU 

•»  2aVi**9*  pairs  are  nuX*red  1-18.',/, 

♦  PM.'Ttw  lirst  laaae  is  nuatur  '.Al,/, 

♦  2CX 'Average  starting  with  pair:  [1]  J,\) 

«A6{*j7S«jMEC  ^  ^ 

m  Fts5t*Tl8H,l«) 

IF<rm£6.l)KKG-l 
IF  K0E6.ls.l.&UCB£G.GT.ll)GOTO  42 
,  MUIE(*.*6} 

848  F0RWT{  ..  ,  .  , 

♦  2IX,'  and  ending  kith  pi)-:  [IS]  ■ ,\) 


32 


siSi, 

IMCcU 

15* 

*  wlr*  It  tt-oltUd.',/. 

■  l* ••••' 

com  t 


,8 

*  r«ttlt  10  prjftis  *,/.  . 

*  ?£$•!  I  *  *s*ws  «**« 

*5  ^  V« ,*««•«««  os*  rr«*  fi) 

Mm&SBIT0 

«* ®  ,M> 

*  11  cce»l«e.*./, 

*  wir 10  e^ll,w4 


*«2  hjvo  ret  b«n  prc«s«d.\/. 

cam  i 

at  mvm 

*  ,M3<  su* (tys,?  u)  **v) 

llikt^MSSS:! 

‘fE3S;;a,S4Sifei]S5F'is'-«»K-'iw) 


aos£  {6ij 

l!UV£.f 

com  \ 

— --c^CcV  staius  fcofare  txlt 
*92  Ftttttt 

*  i?fAj2I5fK?.?Vrrenfc  lW}?  w  '•» 

JFjZ.EiJr^.Z.EQ.yjCOTO  5? 

nrSfflRBl  “ll 

CIOSCJM) 

l«|IE(»,*99) 

*99  FQftn*tC>. 

♦  ?*V|—  ■■  -  -  -  ■  »./ 

4  |S5‘lL  £f°  *  £ISI  WC£  PROCESSOR  (GO) [•!/! 


'  ,rj’  1 1 _  inwvcMyv  VIM/ 

♦  «A.  1,1  1  ■■"  -  .  ■■ 

*iKh 

SOBfOOTlNE  AmG(NO.NOUT,KfiEG,KENO> 

BSSR  SS;kii‘!)':K^jlt,'l!"r 

""Cl'^reglsurj  tod  compute  sw 
w  i  n*l  1 2 

DO  4  J-l.S* 

00  *  1-1.45 
9  fCUT(I.J.H)** 
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CONTINUE  . 

wut« 

CONI  NIC 

RETURN 

ENO 


SUBROUTINE  IGO(N) 
DUUUCTERM  ^ 
REAO(V(AinZ 

lFU.EQ.'C’.OU.EO.Y)?** 

HEtUW 

ENO 


ifftrenew 


SUBROUTINE  JIE5TCIV1 , 1V2.LE5CU.  10CU 

1vi?iv2  data  values  (Input) 

JO^ION  IV2-IV1.  IV1-1V2  dl 
10ft  last  large  ju«i  12SM*J.-J.*! 

I  DEL  sl:e  of  this  Ji»p  (output) 

II  threshold  levels 
—variable  type* 

OURACIER‘1  2 

INIEGf.R*2  IVI,IV2,106,10£L.IDEL.IT(H) 
INTEGERS  JSISn,  |K«.  I  jMP.KSa(2i 
DATA  iT/e.4<.12j.J92.256(j;«,JBt, 

«8.4l2.5?i.MC.7e(,K8.8J2> 
—find  slie  of  Jutp 
K.IV2-1V1 
SIGNM 

f  ioe.tt.8)coio  £ 


GN-l 

ce-ioe 

u«>»>e 

IF(JOC.LT.IT 


.LT.IT 

DC.LT.IT 

C6.LT.IT 


IB 

£8 


14).ANO.IDe.ffi,|T(12))LW*fi 
IZj.WO.lOE.GE.iTdB  IJMP*S 
lflj.ANO.lDC.G6.ITl8  jlJMPM 
8). WO.  IOC, GE. 11(6)  IJMP*5 
10e.LT.IT  6  .WO.I06.G6.ItW  1J«*»2 
IDC.lt. ITd  . WO. IOE,G6.IT(2j  IJMP.l 
I0£.LT.1T(2  .WO. IOC.06. llhiilJf¥»B 
unp.ce.fl.wo.UttP.LE.ncoto 
.  ,  IJNP.G€.3.AN0.1J,TP.IE.5|gOIO  8 
——Illegal  or  unusual  Jumps 
«SG(1)*ISI®«IJ(1P 
KSG  2WDE 
CALL  KSSO(l.NSG) 

iF(Uw.Eg.i4)unP*fl 

IFlUNP.LE.l.«.LDEL.EO.fl)COTO  18 

I0UN*(-1)*ISIGN*L0CL 

I F  ( I0LR1.GT.  8)G0T  □  10 

K&m-isiGNMJnP 

MSG  2  ROCL 

CALL  nCSSG(2,ttSG) 

—return  the  jiatp  value 
IOCL»ISlGN*IJNP 
1 F ( 1 JHP.GT.l )LD€L«10£L 
CONTINUE 
RETURN 
ENO 


SUBROUTINE  MESSG(N.HSG) 

CHARACTER‘1  Z 
INTEG£R*2  «SG(2) 

OATA  Z/'SV 

lF(N.E0.1)WlTE(*,Jfll)«SG 
IFIn.EO.2  WlTEl«,Jfl2  NSG 
301  FCRHATIZSX.SdH"), 

+  '  Jiup  of  ' . 12. *  ( ‘ .14. * >  encounter ed',\) 
382  F0RMAT(25X.5hH*), 

♦  '  Ju«p  ',12,'  follows  jump  ‘,12, \) 

iftzTeq. *6* .6r.z.eo. 'q' iGOTO  § 

REA0(*,'(A1)')Z 
9  RETURN 
ENO 


IWt 


pf* “ 

ws^r 

4  3o,U 

iafe.-- 

Ti«;*tvut»i 
IlMRWtfili 
faTtwx 
JaPM5*T-n« 
Spajaf 


l&nltltllwiH&i'Si 


■nw 


.  its,-  *k«  piraseter*:  . 

*  hr*,‘  win/s*  ,}*,/, 


UpS;,,.,.. 

*r*r*d|w»t  pirmWf* 

}ti  ^n.W(i«,Wrtt  JW4*<UfJ  tV/NI?  (N) 

tl^i’Jrll^X.VJCOio  12 

J14  yo«£{(il*^Ncw  iN-tiKflld:  *,\) 

m  KMiir 

ITWUIOUK 

5,1  ‘A’ 

5,<«  rfVtrrwi  saturation  below  threshold*) 

caio  i 

16  l&MSAT-lll* 


'.M 


ar.istf 

str-sti*, 


- */)2G. 

-—convert  for  printing 
12  00  II  J*l,$fl  r 
00  16  1*1.45 

inw.LT.jita  isv-ii»« 

IFiW.GI.iSATjNV.iSAT 

MMW-IIW 

V‘hV 

V.V/SLI 

‘hV«Y 

JR£G*IT0P-f.7 

IFjjREG.GT.jTOPMREG'lTCf 
IFC IKCG.LT. 180T) IRCC* I80T 
f*(l.J)*REG{l) 

12*  1-1 5*2 
J2*(J*1  j*2 
00  14  Jl*l,2 
JJ»J2*J1 
00  14  11*1,2 

u*i2«n 

!iM'RK(u  „ 

™*H  tnd-of-rtcord  earks 
lR£G*fWi< 


P(46,J)*REG(1) 


17 

18 


JJ*J2vJl 


ill! 


RETURN 
EM! 
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m«  *4te 

WMKM.flW 

fc 

*4  rHaiHjt 

—•«•{«.  «u  t*f)t  (c*l  HJ(  f+l  *■«) 

|*U  U  fdt  tV/M)?  (K)  ',M 

’1* 

**  SMjiulL’rtU 

Til  rcwut(7;>(;K-)J 

M  CSNIJNUE 

««  *tnm 

tNO 

suwoun*  wccKtyj.m.K) 

CWUUCU*M  2 . .  , 


*-•-••— -fC-nt'.Vvt 


RUX-jJTtJ 

MiMii 


*—• — *«Upl»y  im?<  fMitf 

Jli  rcwAtt/.MXj’CfsrfiiUftj  1**5J  '.U)„ 
*........r<cpflUruct  ttAWf  co!w«i  (M*ro*  Hall) 


mi 

a  m  J*i.M 
>u 

!J»J“I 

Vl^NGU.JJ.i:) 

k®j. 


tv* .  iv2.ta 

*  gsilarl^"ht 

. fr  rlaht  Sl& 

00  J*1.U 

NM 


IW.tOCUOCL) 


00  M  J*2MS 
1*1 

] 1 1*11-1 
vi.wimhj.*: 
iv2«fohl,J,Kj 
OIL  Jifciljlvi. 


1V2,[.KI,I0EU 


Wl  JUiitdVi.lVS.im 
jcui.Jwiodil.JHDa 
continue 


33  CONTlfJUE 
«  CCNIINUC, 

... — f*,  left.  s!tf< 
00  M  J*1,M 
KM 
lCEl*« 

icel*« 

00  41  1*1.22 
it-23-: 

111*11*1 


ivj^ojn.j.Ki  _ _ 

call  jlEsUlvi.m.iocuoeu 
w(U.j)*todii.J)*!oeL 

4*  CONTINUE 


5i  CONTINUE 
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m 


10  cartel  iw  o*u 


iO* 


5i 

***•  ati/Diu  f«“*  M<  hsvt  been  comets^.') 
•Mm  mm  »M  oin 


U 


mx«nuD( 

RSdW^tWI  , 

illiila 


«ux*w 

ntfHsi 


!i:il 


ffi  ^  t^W 


jj  *W|;.ss*> 

'  »  ut’lhrf  offset  to  be  Mde*  to  Mta):  \\) 

,  *aa{*,ffit)Tetn 

?<<  ^^l^*|^£2,,M,5w<-cr,,l!a,sai0  55 
nuwvx-nW 

BtN»# 

•  2/x, ' v»ig<$  have  been  offset  by  * ,15,/. 

«  UX.'K «w  «<lewa  sm  elninuo  we  \2l5) 
*•••— —rtPlxt  old  cow  *U*>  new 

*  £fcV* 

CO  MJM.J 
it  CONt 


Nt 

"w”c&<tiNue 

KCUKN 

VC 


MX 
t 


StJWOyilfC  SU8JBIMJ) 

/  :  :: 

Nl*-1 


|XI£C£R*J  NQH5, 


N> 


00  9  JMHM.lt 
KWtMlfM)*2*l 
K2*M 


00 

00 


m 


Wf!.J,K2)*!fl{J,J.K2)*N0(I.J.Kl) 
NOil  |JiKt)»NI»f>0{l,J,K2) 

9  CCNllMX 
KTURM 
DO 
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